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Abstract 

We consider four lepton-number violating (£) processes: (a) neutrinoless double-beta decay 
(Of/3/3), (b) AL = 2 tau decays, (c) AL = 2 rare meson decays and (d) nuclear muon-positron 
conversion. In the absence of exotic Jj interactions, the rates for these processes are determined 
by effective neutrino masses {m)^ 2 , which can be related to the sum of light neutrino masses, the 
neutrino mass-squared differences, the neutrino mixing angles, a Dirac phase and two Majorana 
phases. We sample the experimentally allowed ranges of {m} e ^ based on neutrino oscillation 
experiments as well as cosmological observations, and obtain a stringent upper bound (m)^ 2 <j 
0.14 eV. We then calculate the allowed ranges for (m)^ ^ from the experimental rates of direct 
searches for the above AL = 2 processes. Comparing our calculated rates with the currently or 
soon available data, we find that only the 0^/3/3 experiment may be able to probe {m) ee with a 
sensitivity comparable to the current bound. Muon-positron conversion is next in sensitivity, while 
the limits of direct searches for the other AL = 2 processes are several orders of magnitude weaker 
than the current bounds on (m)g g ■ Any positive signal in those direct searches would indicate 
new contributions to the interactions beyond those from three light Majorana neutrinos. 
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I. INTRODUCTION 



Fermion masses and flavor mixing are among the most mysterious problems of contempo- 
rary particle physics and they have posed major challenges to particle theory and experiment 
for decades. Further understanding of these issues should eventually shed light on fundamen- 
tal phenomena like CP violation, flavor- changing neutral currents, baryon-number (B) and 
lepton-number (L) asymmetry in the Universe and will hopefully lead to a more satisfactory 
unified theory of flavor physics ^J. 

In the Standard Model (SM) of strong and electroweak interactions, neutrinos are strictly 
massless due to the absence of the right-handed chiral states (u R ) and the requirement of 
SU(2)l gauge invariance and renormalizability. Recent neutrino oscillation experiments 
have conclusively shown that neutrinos are massive [2]. This discovery presents a pressing 
need to consider physics beyond the Standard Model. It is straightforward to introduce a 
Dirac mass term m D (T7[i> R +h..c.) for a neutrino by including the right-handed state, just 
like the treatment for all other fermions via the Yukawa couplings to the Higgs doublet in 
the SM. However, a profound question arises: Since u R is a SM gauge singlet, why should 
not there exist a gauge-invariant Majorana mass term \Mv R v R in the theory? In fact, 
there is a strong theoretical motivation for the Majorana mass term to exist since it could 
naturally explain the smallness of the observed neutrino masses via the so-called "see-saw" 
mechanism m u « m 2 D /M 0. From a model-building point of view, there are many scenarios 
that could incorporate the Majorana mass. Examples include R-parity violating interactions 
(AL = 1) in Supersymmetry (SUSY) pj, Left-Right symmetric gauge theories (Hf, grand 
unified theories p, models with exotic Higgs representations @, Hf and theories with extra 
dimensions |9]. One may also consider constructing generic neutrino mass operators to 
parameterize the fundamental physics effects in a model- independent manner [lOj. 

Besides the phenomena of neutrino flavor oscillations and possible new CP-violating 
phases, the Majorana mass term violates lepton number by two units (AL = 2), which 
may result in important consequences in particle physics and cosmology. Although the pre- 
vailing theoretical prejudice prefers Majorana neutrinos, experimentally testing the nature 
of the neutrinos, and lepton-number violation {!£) in general, is of fundamental importance. 
The basic process with AL = 2 is mediated by 



W~W~ 



'1 C 2 5 



where the W~ are virtual SM weak bosons and = e, //, r. By coupling fermion currents 
to the W bosons as depicted in Fig. [H and arranging the initial and final states properly, one 
finds various physical processes that can be experimentally searched for. The best known 
example is the neutrinoless double-beta decay (Of/5/3) [U, 13; which proceeds via the 



parton-level subprocess dd — > uu W~*W~* — > uu e~e~ . Other interesting classes of L 
processes involve tau decays such as r~ — > e + (/i + )7r~7r~ etc. and hyperon decays such 
as T>~ — > S + e~e~, S~ — > p\i~ \T~ etc. |15I |. One could also explore additional processes like 
e~ — > n + conversion . One may also consider searching for signals at accelerator and 
collider experiments via e~e~ -> W~W~ [TJ, e^p -> u e (u e )£f£fX Q, u e (u t )N -> PififX 

Assuming no additional contributions from other exotic particles that have h 1 interactions, 
the matrix element for h 1 processes is proportional to the product of two flavor mixing matrix 
elements and a h 1 mass insertion from a light Majorana neutrino 

= 1X^1*^2^1- 
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FIG. 1: A generic diagram for AL = 2 processes via Majorana neutrino exchange. 

The (fn) e £ are called "effective neutrino masses". Experimental searches for ft processes 
will directly measure the effective neutrino masses squared, and thus probe the fundamental 
parameters of neutrino mixing angles and phases and their masses. 

In this paper we study AL = 2 processes related to W~W~ — > tyt^ ■ We establish our 
conventions and lay out the general expressions for the effective neutrino masses in Sec. |HJ 
We then calculate those quantities, based on the current knowledge from atmospheric, so- 
lar and reactor neutrino oscillation experiments. We include the constraints on neutrino 
masses from a joint analysis of Wilkinson Microwave Anisotropy Probe (WMAP), Cosmic 
Microwave Background (CMB) data, the Sloan Digital Sky Survey (SDSS) large scale galaxy 
survey and bias, the SDSS Lya forest power spectrum and the latest supernovae SNIa sam- 
ple. We thereby limit the allowed ranges of the {Tn) e e . In Sections |HIl IIVI fVl and IVII we 
study the four AL = 2 processes, 

(a) neutrinoless double-beta decay (0v{3/3), 

(b) JJj tau decays r~ -> £ + MfM 2 ~, 

(c) # meson decays Mf -> M^it^t, 

(d) nuclear muon-positron — e + ) conversion. 

We calculate the transition rates and determine the allowed ranges of (m)^ # from the bounds 
set by direct experimental searches for the above AL = 2 processes, and then compare the 
results with the bounds obtained based on the neutrino mixing inputs in Sec. HI1 We find 
that the upper bound from Oz//5/3 is most sensitive to the model-parameters and is at the 
same level as the constraint from Sec. |HI The bounds from the other three classes of 
measurements are significantly weaker although they probe different combinations of the 
parameters. In the future should we observe a ft signal in one of those channels, it would 
indicate non-standard physics beyond the contributions of light Majorana neutrinos. We 
draw our conclusions in Sec. IVII1 Some technical details in calculating the transition rates 
for the above processes are presented in the Appendices. 
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II. EFFECTIVE NEUTRINO MASSES 



In terms of neutrino mass eigenstates Vi the charged current interaction Lagrangian is 
written as 

c cc = —% E E ww.+i-c. (i) 

V 2 £=e,/i,ri=l,2,3 

where Px, is the left-handed projection operator (1 — 75)/2, is the Maki-Nakamura- 



Sakata-Pontecorvo (MNSP) mixing matrix element |2jj between the lepton mass eigenstate 
£ = e, fi,r and the i th neutrino mass eigenstate. It is conveniently parameterized by 



V 



iS i^- i& i5 i^2. is 
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in the notation q = cos#j and s, = sin#j. The mixing angles 9 a and # s are relevant to 
atmospheric and solar oscillations, respectively, and the angle 9 X is presently unknown, 
except that it is bounded by the CHOOZ reactor data. The phase S is a Dirac phase and 
02, 03 are Majorana phases. 

The general subprocess of // is neutrinoless dilepton production from two virtual W 
bosons as depicted in Fig. [T] 

W-W-^WZ, (3) 

which can occur only if neutrinos are Majorana particles. This subprocess changes the 
lepton number from L = to L = 2 and the observation of J/i would establish that neutrino 
is a Majorana particle. The leptonic subprocess of Eq. (j3J) occurs via Majorana neutrino 
exchange and is given by the product of two charged currents 

i 

As presented in Appendix the transition rates for light neutrino exchange are proportional 
to the squares of effective neutrino masses (m)^ defined as 

( m >^ 2 = lE y wV^ m * |, (5) 

i 

where £±,£2 — e,fj,,r, i = 1,2,3. Since the mixing factor is symmetric Vfu Vf^j = V^V^i, 
there are six different (m)^. The explicit expressions for {w) t t are given 22] by 

m) ee = \m 1 clcl + m 2 s 2 s c 2 x e i4 ' 2 +m z s 2 x e i4, % 

m )e U = \mic s c x (-s s c a - s a s x c s e lS ) + m 2 s s c x (c a c s e %<t>2 - s a s s s x e l{<t>2+S) ) + m^SaS^e 1 ^ 5 ^, 



m) eT = \mic s c x (s a s s - s x c a c s e %s ) + m 2 s s c x (-s a c s e t(t ' 2 - s s s x c a e %{:l)2+S] ) + m 3 s x c a c x e 1 ^ 
m ) uu = \mi{s s c a + s a s x c s e iS ) 2 + m 2 {c a c s e i ^ - s a s s s x e l ^ + ^f + m 3 s 2 a c 2 x e l ^ +2 ^\ 



1 • 



m ) M r = \m 1 (s s c a + s a s x c s e l )(s x c a c s e l - s a s s ) + m z s a c a c x e 



2 i(03+2«5) 



+ m 2 (s a s s s x e l ^ 2+5) - c a c s e i4,2 )(s a c s + s s s x c a e iS ) 



1 1 



( m ) = \ mi (s a s s - s x c a c s e iS ) 2 + m 2 {-s a c s e i ^ - s s ^c Q e i( ^+ 5) ) 2 + m 3 c 2 a c 2 x e i{,t>3+25) \. (6) 
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We can see that the {m) eii2 are functions of the oscillation angles (9 a , 9 S and 9 X ) and the 
neutrino masses mj. The three neutrino masses can be expressed in terms of other three 
measured quantities: the sum of neutrino masses and the two mass-squared differences, 

E = mi + m 2 + m 3 — m\-\- \Jm\ + 5m 2 s + Jm 2 + 5m*. (7) 
The atmospheric (a) and solar (s) mass-squared differences are defined as 

5m 2 a = m\ — m\, 

5m 2 = m 2 — ml, (8) 

where 5m 2 > for the normal hierarchy (NH) and 5m 2 < for the inverted hierarchy (IH). 

The above expressions provide a convenient formalism to study the range of values of 
( m )hi2 as f uric ti° ns °f the angles and the phases (both Dirac and Majorana). In particular, 
we study E versus (iTi) iii2 for both the normal and inverted hierarchies. To do this in a 
comprehensive manner, we carry out a Monte Carlo sampling of the oscillation parameters 
and phases. The Dirac phase (5) and the two Majorana phases (02 and (f) 3 ) are allowed to 
range between and 2ir. The atmospheric oscillation data gives bounds on [23] 

1.9 x KT 3 eV 2 < \5m 2 a \ < 3.0 x 10~ 3 eV 2 , sin 2 29 a > 0.9 at 90% CL. 

In fact, the bounds on 9 a vary with 5m 2 a and in our computation these bounds are obtained 
from the 90% CL 5m 2 , versus sin 2 26 l a plot of Ref. [24|, which was obtained in an L/E analysis 
of only selected high resolution FC (fully contained) and PC (partially contained) events. 
The analysis of the full data set from the same running period |25| gives slightly different 
constraints; the constraints on sin 2 26 l a are slightly better from the full data set but the L/E 
analysis better constrains 5m 2 . Similarly, for the above range of 5m 2 , reactor data places 
bounds on 

sin 2 ^. < 0.06, 

which also vary with 5m 2 . We use the limits obtained from the CHOOZ 90% CL exclusion 
plot of Ref. [26]. Finally a joint analysis of the solar and reactor oscillation data limits the 
parameter ranges to 

7.6 x 10~ 5 eV 2 < 5m 2 < 9.1 x 10" 5 eV 2 , and 0.31 < tan 2 9 S < 0.52 at 90% CL. 

The bounds on 9* yary with 5m 2 and these are obtained from the 90% CL 5m 2 versus 
tan 2 ^ s plot of Ref. [27]. All the inputs to the Monte Carlo sampling are summarized in 
Table I. We can further constrain the range of (fn) e e by imposing limits on S obtained 
from cosmology. The current best limit 

E < 0.42 eV at 95% CL (9) 

was obtained from an analysis of WMAP, the SDSS galaxy spectrum and its bias, the SDSS 
Lya forest power spectrum and the latest supernovae SNIa sample, assuming a spatially 
flat universe and adiabatic initial conditions [28j. Throughout this paper, we will adopt the 
cosmological bound of Eq. (jUJ). We would like to point out that more conservative analyses 
without including the Lya forest power spectrum exist. These lead to larger values of E, 
such as E < 0.75 eV @ and E < 0.54 eV at the 2a level. 

The results of the Monte Carlo sampling along with constraints from cosmology are shown 
in Fig. |21 for both normal and inverted hierachies by solid and dashed curves respectively. 
We first point out the following general features: 
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TABLE I: Inputs to Monte Carlo sampling for calculating effective neutrino masses 



Parameter 


Input 


\5m 2 a \ 


1.9 x 10" 3 eV 2 - 3.0 x 1(T 3 eV 2 [23] 


5m 2 s 


90% CL 5m 2 s versus tan 2 6» s plot [27] 




90% CL 5m 2 a versus sin 2 20 a plot [24] 


e s 


90% CL 5m 2 s versus tan 2 6» s plot [27J 


Ox 


CHOOZ 90% CL exclusion plot [2£] 


5 


to 2vr 




to 2vr 


4>3 


to 2tt 



For large values of the minimum neutrino mass, typically m min > 0.1 eV, the mass 
differences are unimportant and thus m 1 ps m 2 ~ m 3 or S > 3m min . Due to unitarity 
of the mixing matrix, we can see that the maximum values of effective masses obey 
(m) ei£2 « m m in %j S/3. On the other hand, for smaller values of m m i n , S is governed 



by the larger mass difference, namely £ > y 5m? a ps m 3 « 5 x 10 2 eV for the NH 
scenario, and £ > 2 J bm\ « 2m! « 9 x 10~ 2 eV for the IH scenario. 



• The cosmological observations of £ < 0.42 eV put an upper limit on (m)^ at about 
0.14 eV. 

• The normal and inverted hierarchies are indistinguishable at the current level of sen- 
sitivity. When E is smaller than the present cosmology limit, the difference (in upper 
limits) between normal and inverted hierarchies may become significant as is apparent 
from all the plots. We are just above the region where the NH and IH results begin to 
separate. An improvement in the accuracy of either of the two observables by a factor 
of two or better would begin to provide a sensitive probe to distinguish the NH and 
IH scenarios. 

We note the qualitative difference between Fig. Ufa) and the other panels in Fig. |21 The 
allowed region for (m) is between the curves for both normal hierarchy (solid lines) and 
inverted hierarchy (dashed curves). It is more stringent than the other (m)^ 2 because (fn) ee 
depends on two fewer parameters (9 a , 5) than the others. Also the specific combination of 
oscillation parameters for (m) does not lead to complete cancellations or vanishingly small 
contributions (unless for a small range of E in the NH scenario) unlike others. Our results 



for the range of values of {m) ee are similar to the analyses previously done in Refs. [31, 
and an updated version presented in Ref. 13311 . Similar analyses for the other {m) e e for 
specific scenarios were considered in Ref. [34]. We will not show the results for (m) and 
(m) TT since they are qualitatively very similar to (m) er and (m) , respectively. 

By setting limits on the f, process decay rates and cross sections, direct experimental 
upper bounds on effective neutrino masses can be obtained as to be discussed in the next 
sections. We then compare them with the results from this section. 
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FIG. 2: (a) Upper left: allowed regions between the curves for (m) ee versus E; (b) upper right: 
allowed regions to the right of the curves for {m) e ^ or {m) eT ; (c) lower left: same as (b) but for 
( m ) fin or ( m ) T r' (°0 lower right: same as (b) but for (m) . The bound for £ at 95% CL from 
cosmology is shown in all figures by the vertical line. 



III. NEUTRINOLESS DOUBLE-BETA DECAY (0i//?/3) 

The decay rate for neutrinoless double-beta decay (QufiP) is proportional to (m) ee ; (m) ee 
is plotted in Fig. |2fa). The theoretical formalism for calculating the decay rate is given 
in Appendix El Ref. ^2 summarized the latest experimental limits of Ov/3/3 for various 
isotopes; the results are reproduced in Table El The experimental bound on (m) ee has 
improved from 5 eV in 1992 to about 1 eV in the most recent experiments. The best limits 
come from the two 76 Ge experiments which are Heidelberg-Moscow and IGEX respectively. 
Although 0^/3/3 is a vital experiment to determine the Majorana nature of neutrinos, we 
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TABLE II: Experimental bounds on half-life time of for various isotopes from Ref. [12j and 

the implied upper bounds on {m) ee . The last row lists the cosmological bound as obtained in the 
previous section. 



Isotope 


Hali-liie (years) 


Hee ( eV ) 


Year of published paper 


48 Ca 


-i A -I r\99 

> 1.4 x 10 2Z 


< 7.2 - 44.7 


2004 


(b Ge 


> 1.9 x 10 25 


< 0.35 


2001 


7fi r-\ 

' b Ge 


> 1.6 x 10 Zi> 


< 0.33 - 1.35 


2002 


' D Ge 


= 1.2 x 10 25 


= 0.44 


2004 


82 Se 


> 2.7 x 10 22 


< 5 


1992 


100 Mo 


> 5.5 x 10 22 


< 2.1 


2001 


116 Cd 


> 1.7 x 10 23 


< 1.7 


2003 


128 Xe 


> 7.7 x 10 24 


< 1.1 - 1.5 


1993 


iso Te 


> 5.5 x 10 23 


< 0.37-1.9 


2004 


136 Xe 


> 4.4 x 10 23 


< 1.8 - 5.2 


1998 


i50 Nd 


> 1.2 x 10 21 


< 3 


1997 


Cosmology 


none 


< 0.14 


this paper 



note that the uncertainty in nuclear matrix elements would result in an uncertainty as large 
as a factor of three in the inferred value of (m) ee from an observation of the decay process 

Hi. 



We list our result in the last row of Table ITDbased on Fig.|2fa) which is largely determined 
by the cosmological bound. We see that our result is slightly stronger than the current 
experimental bounds. The bounds we obtain for {m) ee are based on the Monte Carlo scan 
over the fundamental parameters in the neutrino sector as given in Table H] and thus are 
independent of the nuclear matrix elements. The large uncertainty comes in when we predict 
the decay rates for the nuclear isotopes. 

A recent publication claims evidence for 0V/3/3 at the 4.2a level |35j], but the result is 
controversial. In Ref. |35J], {m) ee was determined to be in a range between 0.2 to 0.6 eV at 
99.73% CL for a particular choice of the nuclear matrix element, and becomes 0.1 to 0.9 eV 
if allowing a ±50% uncertainty of the nuclear matrix element. With the bound on £ at 95% 
CL, the first range is disfavored by the limits we obtained and the second range allowing a 
larger uncertainty leaves a very narrow range for {m) ee . 

Tritium /3 decay experiments also probe the absolute scale of the neutrino mass. The 
current limits from cosmology are better by an order of magnitude compared to the tritium {3 
decay limits |33| . Future limits from the KATRIN tritium /3-decay experiment are expected 
to be 0.30 eV (0.35 eV) at the 3a (5a) level. The present limit on (m) ee from cosmology is 
still stronger than this expected improvement from KATRIN, but the KATRIN experiment 
will provide an important direct confirmation. 
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IV. LEPTON-NUMBER VIOLATING TAU DECAY 



In this section we examine tau decays into an anti-lepton and two mesons 

r-^tM^M^, (10) 

which is a process with AL = —2. The relevant effective masses are {m) eT and (m) as 
shown in Fig. 12(b) and Fig. |2Jd) respectively, from the current parameters from neutrino 
oscillation experiments. The constraint from cosmology of Eq. © gives an upper limit 
of 0.14 eV for (m) eT and (m) . In Appendix we give the calculations for the decay 
branching fraction of the process (|10|) in terms of (fn) 2 £T . We express the branching fraction 
in an intuitive form as 

m- 33 \ V ckm v ckm i2 f fat! Im 2 V / 1777 MeV \ 2 / (m) e V 

where $ is the phase space integral over the squared matrix element and can be evaluated 
numerically. For small values of (m) the branching fraction induced by a light Major ana 
neutrino is seen to be very small. 

A direct search for neutrinoless tau decays has been made at the CLEO II detector 
at Cornell Electron Storage Ring (CESR). Twenty eight different decay modes have been 
studied and the limits on the branching fractions were reported in [36|. The experimental 
limits for various decay modes are typically of the order of 10~ 6 , as given in Table 11111 
From those, one can determine upper bounds on (m) er from Eq. (jllj) . as given in details in 
Appendix Unfortunately, the obtained bounds are very weak compared to our inferred 
cosmology bounds less than 1 eV, as shown in the last column in Table IIHI In fact, the 
current formalism for calculations in terms of the effective neutrino masses is valid only 
for the light Majorana neutrino exchange when the mass is much less than the energies 
available in the reaction. Thus the entries with such large values in this Table lose the 
original meaning of the effective neutrino mass. We nevertheless include these values here 
and henceforth to indicate how much improvement would be needed to be sensitive to the 
light neutrino contributions. This information is useful to see which process may be more 
sensitive to what operator and to what extent. On the other hand, any observation of a 
I/j signal in these channels at the current values would strongly imply contributions beyond 
those of light Majorana neutrinos. Hence it is important not to neglect the experimental 
study of these processes even though the limits seem too weak in comparison with the 
Majorana neutrino mechanism. 



V. RARE MESON DECAYS 

We now investigate the J/j processes in which a meson decays ji^, Ht| into another meson 
and two like-sign leptons 

M+ - M 2 "£+4. (12) 

These decays are similar to the tau decay modes described in the last section. For the various 
decay modes, the effective neutrino masses involved are {m) ee , (m) and (m) depending 
on the final state leptons. Again, we plot their allowed values based on the known neutrino 
parameters, as shown in Fig. 12(a) for {m) ee , in Fig. |2Jb) for (m) (indistinguishable from 
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TABLE III: Experimental bounds on branching fractions in AL = 2 tau decays from 36] and the 
implied sensitivity to probe the corresponding effective neutrino masses. 



Decay mode 


Bexp 


(m) er (TeV) 


r~ — * e + ir~ir~ 


1.9 x i(r 6 


12 


t~ — * e + ir~K~ 


2.1 x 10~ 6 


46 


T~ -> e+K-R- 


3.8 x 1(T 6 


730 


t~ —* /i + 7r~7r~ 


3.4 x 1(T 6 


20 


r~ —* fi + ir~K~ 


7.0 x 10~ 6 


100 


t' -» ij+K-K- 


6.0 x 10~ 6 


1000 



(m) er ), and in Fig. |2fc) for (m) (indistinguishable from (m) Tr ). We infer a generic upper 
limit for (m) l ^ to be 0.14 eV from constraints from cosmology Eq. (jHJ). The branching 
fraction for the rare meson decay modes is 

DP ~ m-29 \jrCKMjrCKM\2 ( T Mi \( fkhfhh \ ( m M x \ 3 ( ( Tn )l 1 i 2 \ ~, nq \ 

BR ~ 10 ^ V * 1 U.OxlO-8 s) {(100 MeVy) {iGeVj {Te^ ) $ ' (13) 

where $' is the phase space integral over the squared matrix element and can be evaluated 
numerically. 

Searches for rare meson decay modes have been made in numerous experiments. Table ITVl 
summarizes the current experimental limits on branching fractions given by . From these, 
direct search limits can be determined on effective neutrino masses and some associated 
calculational details are given in Appendix [Dl Again, the bounds obtained are still much 
weaker than the cosmology bound. Although the K + decays yield the most sensitive bounds, 
they are still many orders of magnitude away. We include the obtained values in Table II VI 
to indicate how much improvement would be needed to be sensitive to the light neutrino 
contributions. There are no direct search limits obtained for (m) TT from the processes 
discussed. However only very weak constraints for BR(B — > Xr + r~) < (9(5%) exist in a 
theoretical analysis [39|. The similar signature B + — > M~t + t + is a possible decay mode 
that would bound (m) TT and should be pursued, but any such bound will not be competitive 
with the cosmology limit unless there is a contribution from new physics beyond the light 
Majorana neutrinos. 

VI. MUON POSITRON CONVERSION 

The nuclear muon to positron conversion process is another AL = 2 process that is very 
similar to 0vf3f3. When a muon propagates through matter, ordinarily it interacts with a 
proton in a nucleus and produces a neutron and a neutrino, which is similar to inverse beta 
decay. However, if the neutrino is a Majorana particle, it is possible that a muon can interact 
with two protons and produce two neutrons and a positron. The leptonic part of the decay 
amplitude is exactly the same as tau decay and the nuclear part will lead to nuclear matrix 
elements analogous to 0v(3f3. The fundamental interaction is parameterized by (m) and the 
current bound from neutrino oscillation experiments and cosmology is shown in Fig. E[b). 
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TABLE IV: Experimental bounds on branching fractions in AL = 2 rare meson decays [38j and 
the implied sensitivity to probe the corresponding effective neutrino masses. 



Decay 


mode 


Bexp 


{m) Mo (TeV) 


K+ -» 




6.4 x 10" 


-10 


0.11 


K+ -» 


7T~/i + /i+ 


3.0 x 10" 


-9 


0.48 


K+ -» 




5.0 x 10" 


-10 


0.09 


,D+ 


7r~e + e + 


9.6 x 10" 


-5 


320 






4.8 x 10- 


-6 


76 


,D+ -» 




5.0 x 10" 


-5 


170 


L>+ -» 




1.2 x 10" 


-4 


1900 


L>+ -» 




1.3 x 10- 


-5 


670 






1.3 x 10" 


-4 


1500 






6.9 x 10" 


-4 


200 






2.9 x 10- 


-5 


42 




7r~e + /i + 


7.3 x 10" 


-4 


150 




K-e+e+ 


6.3 x 10" 


-4 


990 


^, + - 




1.3 x 10- 


-5 


150 




K~e + fj, + 


6.8 x 10" 


-4 


740 


£?+ -► 


n~e + e + 
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-6 


420 
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1.4 x 10- 


-6 


400 


B+ -» 


7r~e + /i + 


1.3 x 10" 


-6 


270 


5+ -» 




1.0 x 10" 


-6 


1300 


-> 




1.8 x 10- 


-6 


1800 


5+ -> 


K~e + n + 


2.0 x 10" 


-6 


1300 



An experimental bound on the branching ratio of muon to positron conversion on titanium 
was reported in j40j 

B = r <" + ; ^ e+ ±%> < 1.7 x 10-". (14) 
T(Ti + fji — > + be) 



The experimental limit on (m) was obtained from this branching ratio limit by Ref. [41 
to be 

(m) e/J < 17 (82) MeV, (15) 

where the created proton pairs are in singlet (triplet) state. Although still larger than the 
bound from oscillation plus cosmology, this can be the next most sensitive probe to the 
Ifj processes after 0v(3/3. However, others 42] argue that the theoretical expression for the 
decay rate was overestimated and they obtain a much lower branching ratio prediction, 

E = 1.60xlO- 25 (^M , (16) 

\ Tflp J 



n 



which would only lead to a weak bound of (m) < 1.3 TeV. Thus there is a large disparity of 



~ 10 in the literature about the inferred limits on (m) from the muon-positron conversion 
process, due to the different treatments of the nuclear transition matrix elements. Beside 
this difference, there is a large uncertainty in muon-positron conversion due to the effects 
from nuclear physics. As a competing channel, the limit from K + decay of 90 GeV is more 
constraining than this 1.3 TeV limit but weaker than the optimistic result above. Moreover, 
the hadronic matrix element for the kaon decay should be better known than that involving 
nuclear physics. 

Another process similar to muon-positron conversion that has not been studied exper- 
imentally so far is nuclear muon capture: /i~ + (Z,A) — ► /i + + (Z — 2, A). This process 
was first studied theoretically for 4A Ti and a branching ratio of 5.0 x 10~ 24 was obtained by 
considering an effective neutrino mass (m) of 250 keV j43[. By including our limit from 
cosmology of (m) < 0.14 eV, we can deduce a branching ratio < 1.6 x 10 -36 . Ref. [44| 
claims that the i mag inary part of the nuclear matrix elements which plays a dominant role 



was neglected in |43( which led to an overestimation. They obtain a branching ratio 



D 1.0 x l0- 23 (^^ . (17) 



With our limit of {m) < 0.14 eV this translates to a branching ratio < 0.75 x 10 . The 



most recent paper on this topic |45J] claims the branching ratio is 10 -6 smaller than the one 
estimated in |43j and would lead to a branching ratio < 10~ 42 if we consider (m) < 0.14 
eV. Evidently there is a disagreement in the literature about the limit on the branching 
ratio by 6 orders of magnitude. Only an effort similar to that for Qvvfi can improve the 
situation. The /i~ — > /i + conversion process will be studied experimentally at the PRISM 
facility and is expected to achieve sensitivity of ~ 10 -13 to this process on 44 Ti nucleus after 
one nominal year run (46[. But this is much below the predicted branching ratios even for 
the most optimistic scenario and will not be accessible in the near future if the process is 
mediated by light Majorana neutrinos only. 



VII. CONCLUSION 

The observation of a process would show that neutrino is a Majorana particle. In the 
absence of exotic Jj interactions, the rates for these processes are determined by effective 
neutrino masses (m)g g , as functions of light Majorana neutrino masses and the mixing 
parameters. We first sampled the experimentally allowed ranges of (Tn) e ^ based on the data 
from neutrino oscillation experiments as well as cosmological observations, and obtained a 
stringent upper bound {m) e e ~ S/3 < 0.14 eV. This cosmology limit is expected to improve 
with a future sensitivity down to X m 0.1 eV [47]. As the limits on £ improve new bounds 
on (iTi) iii2 can be deduced as seen from our plots. In particular, the normal hierarchy and 
inverted hierarchy scenarios may be experimentally differentiated. 

We considered four lepton-number violating processes: (a) neutrinoless double-beta decay 
(0uj3j3), (b) AL = 2 tau decays, (c) AL = 2 rare meson decays and (d) nuclear muon- 
positron conversion. After evaluating the transition rates for these processes, we translated 
the current experimental bounds from direct searches into limits on (m)^ 2 . The best limits 
obtained from experiments on these processes were compared with the cosmology limits 
in Table |V] The Ov/3/3 process is the only process which can currently give interesting 
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TABLE V: Summary of experimental bounds and cosmology limits on effective neutrino mass. The 
lowest bounds on each component of effective neutrino mass are presented with the corresponding 
experiments for these bounds. 





Cosmo bounds on (m)^^ 


Exp bounds on (m)^^ 


Corresponding experiments 


ee 


0.14 eV 


0.33 eV 


Ov00 


e/x 


0.14 eV 


17 MeV (90 GeV)t 


ji~ — e + conversion 


er 


0.14 eV 


12 TeV 


T~ — ► e + 7T~7r~ 


nn 


0.14 eV 


480 GeV 


K+ -► ir-fi + n + 


fiT 


0.14 eV 


19 TeV 


T~ — ► jJ, + TT^7T^ 


TT 


0.14 eV 


none 


none 



t The conservative limit comes from X + — > it e + /x + which, unlike ti — e + conversion, 
does not involve the large uncertainties from nuclear matrix element calculations. 



experimental limits on (m) e e . We note that while experimental limits for Ov/3/3 involve 
large theoretical uncertainties from nuclear matrix element calculations, our cosmology limit 
is independent of any such uncertainties. The other processes have very weak experimental 
limits, that essentially do not impose any meaningful bounds on (Tn) e e . The entries in the 
Tables are only meant to suggest the level of improvement needed in order to sensitively 
probe the light Majorana neutrino mass. On the other hand, the predicted small rates could 
provide a window of opportunity for observing exciting new physics. Any positive signal in 
those direct searches would indicate new contributions to the $ interactions beyond those 
from the three light Majorana neutrinos. 
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APPENDIX A: GENERAL AMPLITUDE OF AL = 2 PROCESSES 

The charged current interaction lagrangian in terms of neutrino mass states is 

£w = — n= E Vj-fPLViWn + h.c. (Al) 

V * e=e,n,r i=l,2,3 
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where Pl = ^(1 — 75). The leptonic AL = 2 subprocess WW — > £ x £ 2 is induced by the 
product of two charged currents 

M%p « E^^^iT^^^T^^)- (A2) 

i 

which can be rewritten using charge conjugation as 

M? e u p oc Yy^VUhl»PL^l u P R t 2 ). (A3) 

i 

The Majorana neutrino fields can be contracted to form a neutrino propagator, and the 
transition matrix element is thus given by 

= 9 4E v ^Vi^P L )^^(YP R i c 2 ), (A4) 
/ j q m>i 

where q is the momentum exchange carried by the neutrino. The q[ term vanishes due to the 
chirality flip. Including the crossed diagram {JL\ <-> £ 2 ) the leptonic amplitude then becomes 

= ^y^-^-^MYY + 7Y)^2. (A5) 
z j q 

If we only consider the contributions from light Majorana neutrinos, namely q 2 ^> m 2 , then 

= ^I^yy + YY)PrV2 Yy^rrii, (A6) 
2q 2 Y 

and is thus governed by the "effective neutrino mass" 



APPENDIX B: NEUTRINOLESS DOUBLE-BETA DECAY (0v(3(3) 

The decay amplitude for neutrinoless double-beta decay (0u(3j3) can be separated into 
leptonic and nuclear parts, 

iM = (M lep ), u (M nuc r. (Bl) 

The leptonic amplitude is given by (jA6|) . In the non-relativistic approximation for the 
nucleons, the nuclear amplitude evaluated for initial ground state to final ground state 
transitions turns into a sum of Gamow- Teller and Fermi nuclear matrix elements defined as, 

2 2 
M nuc = M GT - %M F = (f\ ]T H(r jk , E)T}rl{a 3 ■ a k (B2) 

9a j,k 9a 

where (/| and \i) are the final and initial nuclear states, g a and g v are weak axial- vector 
and vector coupling constants and the function H called the "neutrino potential" has an 
approximate form given in Ref 12]. The decay rate for 0^/3/3 can be expressed as 

[T1]- 1 = G(AE, Z)\M nud \ 2 (m) 2 ee , (B3) 

where G(AE, Z) is the phase space integral. 

For a detailed discussion, in particular the uncertainties associated with the nuclear 
matrix elements, see [l^, EH . 
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APPENDIX C: LEPTON-NUMBER VIOLATING TAU DECAY 



This mode is cleaner in principle than 0i/(3(3 since the hadronic part does not involve com- 
plicated nuclear structure. For the tree level amplitude, the hadronic part can be expressed 
in terms of the decay constants of the mesons in a model independent way. The box diagram 
includes hadronic matrix elements which cannot be simplified in terms of decay constants 
and needs to be evaluated in a model dependent way. We expect the tree level amplitude 
to dominate and do not include the box diagram. It has been argued that in certain cases 
for rare meson decays sub-leading contributions may be appreciable (20I l48l | . Even in such 
a scenario the difference will not be important at the current level of sensitivities and we 
include the more conservative limit from tree level diagrams only. The tau decays and the 
rare meson decays are crossed versions of each other, hence the above arguments are true 
for both. 

The leptonic part of the subprocess t~ — > £ + W~*W~* is obtained by crossing the ampli- 
tude of W + W + -> if it in (JMD 

M% = ^Mer^^. (Ci) 

z i 9 

Combining the hadronic and leptonic parts, the decay amplitude for 

is given by 

iM = (Muv^MfaM^ + (Mi <-> M 2 ) 

= 2GlV^V^f M J M2 £ VrWumMj^r, + J^)Pnv e ] , (C2) 



where V is the quark flavor-mixing matrix elements for the mesons, are meson 
decay constants. The decay rate is then given by 

r = (1 - ^M^^^v^yM^flJ^M Hi ^ (C3) 

where $ is the phase space integration over the matrix elements squared 

$ = — - / F(pi,qj) dx e dx M (C4) 
mjj. J 

A B 

(pi - iir [Pi - Qi) 1 (pi - Q2) 2 

A(pi, qj) = 8(pi • qi)(p 2 ■ q 2 )(qi • 92) - ^ m lh(Pi • 92) (P2 • 9a) 

- 4m M 2 (Pi • 9i)(P2 ■ 9i) + 2m 2 Ml m 2 M2 (pi ■ p 2 ), (C6) 
B(pi, <&) = 4(p! • p 2 )(9i • 9 2 ) 2 - Afo, 9,)- (C7) 

The integration variables xe and a; A/i are the energies scaled by the mass of the decay particle 

x % = ^, (C8) 
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as introduced in the text book |49|. Numerically, $ ~ 1.641, 0.7787 and 0.1455 for the 
modes 7r7r, ttK and KK respectively, neglecting the mass of the final state lepton. 

In the limit that the final state particles are massless, then the phase space $ can be 
written in a simple form as 

$ = \ I [A(x e - 1) + X ^ - 1} -] dx e dx Ml , (C9) 

m 2 . J (x Ml - l){x e + x Ml - 1) 

where the integration limits are given by < Xt < 1 and (1 — X() < xmi < 1- Note that $ 
presents a mass singularity when all the final state particles are considered massless. 

Normalized to the r decay width T T = G 2 F m^/l92ir 3 , the corresponding branching fraction 

is: 

BR=(1- bM^f-GllV^V^fkfk^ Ml * (CIO) 



- 33 '«™i 2 fl^) 2 (™) 2 (!%'* 



~ m- 14 \\rCKMirCKM\2 ( f Mi /m 2 \ ( ( m ) £ T \ , T) 

~ 10 |Kmi Vm > 1 1,(100 MeVy) 

The meson decay constants, CKM matrix elements and r mass are taken from the Particle 
Data Group (PDG) 

U = 130.7 MeV, f K = 159.8 MeV, \V ud \ = 0.9738, \V US \ = 0.2200. 

APPENDIX D: RARE MESON DECAY 

The rare meson decays 

M+(q 1 )^£+{ Pl ) £ + (p 2 )M^(q 2 ) 
have the same Feynman diagrams as tau decay. The decay amplitude is given by 

iM = 2GlV™ M V™ M f Ml f M2 £ V eil V e2l m t u ei ( ( M + ^ A \ 2 )P R ve 2 } . (Dl) 

l i [oi-pir \q\-p2Y J 

The decay rate is then given by: 

r = (1 - ^^ 2 )^^l^7 M ^f M r/^/^r< (m)l h (D2) 
where $' is the phase space integration with the dimensionless integration variables xg 1 and 

$' = j F'(pi,qj) dx ei dxM 2 , (D3) 



Mi 



A B 
F'iPi, Qj) = 7 u + 7 ^7 N2 + (Pi ~ P 2 )' ( D4 ) 
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where A(pi,qj) and B(pj, q$) are given in Eqs. (|C6j) and (|C7[) . and % 2 are the energies 
scaled by the mass of the decay particle and are given by Xi = 2E i /rriM 1 [49]. To have a 
numerical estimate consider the case when the final state particles are massless. Then the 
phase space $' can be written in a simple form as 

$' = / 4(1 - XM 2 )dx ei dx M2 ~ 0.6667, (D5) 

m Ml j 

where the integration limits are < xi x < 1 and (1 — x^) < xa/ 2 < 1- It is interesting to 
note that the integration $' is finite even in the massless limit for the final state particles, 
unlike the case for $ in r decay, due to the anti-symmetric property of the matrix element 
for the two fermions in the final state. 
The branching ratio is then given by 

BR = r M T = (1 - ^^.J^r^G^I^^I 2 /^/^^ (m) 2 ^ d>' (D6) 
« 10- 29 W™ M V™ M \ 2 ( : , ™\ 5 X 



1.0 x 10" 8 s 



^(lOOMeV) 2 ^ UGeVy \ 1 eV J ' { } 

We have used the following constants from the PDG (38[, the decay constant for B from 
(Hoj l and for D from (H]], for obtaining (wi)^ from the branching fractions for the various 
decay modes 

f D = 202 MeV, f Ds = 266 MeV, f B = 190 MeV; 
\V ub \ = 0.00367, |V cd | = 0.224, (V^l = 0.996; 

r K = 1.2384 x 10 -8 s, r D = 1.040 x 1(T 12 s, r Da = 4.9 x 1(T 13 s, r B = 1.671 x 1(T 12 s; 
m K = 493.7 MeV, m D = 1869 MeV, m Ds = 1968 MeV, m B = 5279 MeV. 
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